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Abstract 

Mittag-Leffler analysis is an infinite dimensional analysis with re¬ 
spect to non-Gaussian measures of Mittag-Leffler type which general¬ 
izes the powerful theory of Gaussian analysis and in particular white 
noise analysis. In this paper we further develop the Mittag-Leffler 
analysis by characterizing the convergent sequences in the distribution 
space. Moreover we provide an approximation of Donsker’s delta by 
square integrable functions. Then we apply the structures and tech¬ 
niques from Mittag-Leffler analysis in order to show that a Green’s 
function to the time-fractional heat equation can be constructed using 
generalized grey Brownian motion (ggBm) by extending the fractional 
Feynman-Kac formula |Sch92) . Moreover we analyse ggBm, show its 
differentiability in a distributional sense and the existence of corre¬ 
sponding local times. 
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1 Introduction 


The Mittag-Leffler analysis is part of the field of research which tries to trans¬ 
fer the concepts and results known from Gaussian analysis to a non-Gaussian 
setting. First steps in this direction were made in |Ito88] using the Poisson 
measure. This approach was generalized with the help of biorthogonal sys¬ 
tems, called Appell systems, see e.g. |Dal91l IADKS96| IKSWY98] . which are 
suitable for a wide class of measures including the Gaussian measure and the 
Poisson measure |KSS97] . Two properties of the measure jj, are essential; An 
analyticity condition of its Laplace transform and a non-degeneracy or posi¬ 
tivity condition (see also |KK99j L In [KSWY98] a test function space {MY 
and a distribution space (A/")”^ is constructed and the distributions are char¬ 
acterized using an integral transform and spaces of holomorphic functions 
on locally convex spaces. The corresponding results in Gaussian analysis, 
see e.g. |KLP+96( IKLS96( IGKS99] , paved the way for applying successfully 
the abstract theory to numerous concrete problems such as intersection local 
times for Brownian motion |dFHSW97] as well as for fractional Brownian 
motion |DQDS08l lOSSll] and stochastic partial differential equations, see 
e.g. |GKU991 IGKSOn] or the monograph |H0UZ96] . 


In Mittag-Leffler analysis, the measures /i^, 0 < /3 < 1, are defined via 
their characteristic function given by Mittag-Leffler functions. It is shown in 
|GJRdS15] that the Mittag-Leffler measures belong to the class of measures 
for which Appell systems exist. Furthermore the (weak) integrable functions 
with values in the distribution space are characterized and a distribu¬ 

tion is constructed which is a generalization of Donsker’s delta in Gaussian 
analysis. 

We further develop the theory of Mittag-Leffler analysis in this work by 
characterizing the convergent sequences in the distribution space This 

is done in Section [2l where we first repeat the main steps for the construc¬ 
tion of a Mittag-Leffler analysis from |G,lRdS15] and we introduce the test 
function space {MY and the distribution space {M)~^. In this way we give 
an approximation of Donsker’s delta by a sequence of Bochner integrals in 
the Hilbert space L^(/i^). 

In Section [3l we demonstrate that generalized grey Brownian motion 
(short ggBm) 0 < a < 2, can be defined in the setting of Mittag-Leffler 
analysis using the probability space (iS'(M), H,/i/j), where iS'(M) are the tem¬ 
pered distributions, B denotes the cylinder a—algebra and /r^, 0 < /3 < 1, 
are the Mittag-Leffler measures. This stochastic process was analysed for 
example in |Sch9nl IMMn9] in a slightly different setting. In Section 0] we 
establish the connection of ggBm to the fractional heat equation. This equa¬ 
tion has already been studied for example in |SW891 IMai951 IKochOl lEKOd] . 
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In |OB09] the solution to the time-fractional heat equation is interpreted as 
the density of the composition of certain stochastic processes. We obtain a 
time-fractional heat kernel by extending the fractional Feynman-Kac formula 
|Sch92] to localized (Dirac delta) initial values and justify in this way that 
ggBm plays the same role for the time-fractional heat equation as Brownian 
motion does for the usual heat equation. In Section [S] we apply the char¬ 
acterizations from Mittag-LefHer analysis to show differentiability of in 
and we prove the existence of ggBm local times using Donsker’s delta, 
see also |dSE15] . 

Summarizing, the core results of this article are: 

• Giving a characterization of convergent sequences in the dis¬ 
tribution space see Theorem 12.121 

• Providing an approximation of Donsker’s delta by a sequence of Bochner 
integrals, see Theorem 12.141 

• Constructing the Green’s function to time-fractional heat equations by 
establishing a time-fractional Feynman-Kac formula, see Theorem 14.31 
and Remark 14.51 

• Showing differentiability of generalized grey Brownian motion in a dis¬ 
tributional sense, see Theorem ESI 

• Constructing local times of grey Brownian motions, see Theorem 15.31 


2 Mittag-Leffler Analysis 

2.1 Prerequisites 

In this section, we repeat the construction of Mittag-Leffler measures as 
probability measures on a conuclear space M' from |Sch92) . First we need 
to collect some facts about nuclear triples used in this paper. For details see 
e.g. [SdlTlllRST^ . 

Let be a real separable Hilbert space with inner product (•, •) and 
corresponding norm | • |. Let AA be a nuclear space which is continuously and 
densely embedded in 'H and let M' be its dual space. The canonical dual 
pairing between M' and M is denoted by (•, •), and by identifying 'H with its 
dual space via the Riesz isomorphism we get the inclusions M <Z 'H (Z M'. 
In particular (/,0 = (/,0 for / G ^ G U. 

We assume that M can be represented by a countable family of Hilbert 
spaces as follows: For each p G N let "Hp be a real separable Hilbert space 
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with norm I'l^ such that M C "Hp+i C Tip C 'H continuously and the inclusion 
'Hpj^i C T-Lp is a Hilbert Schmidt operator. It is no loss of generality to assume 
Mp < Mp+i ^p+i "Ho = "H, I'lg = I'I- The space N is assumed to be the 
projective limit of the spaces ('Hp)pGN, that is TV = flpeN^p topology 

on M is the coarsest locally convex topology such that all inclusions M C "Hp 
are continuous. 

This also gives a representation of M' in terms of an inductive limit; Let 
1-L-p be the dual space of 1-ip with respect to li and let the dual pairing 
between H-p and Up be denoted by (•, •) as well. Then H-p is a Hilbert 
space and we denote its norm by |•|_p. It follows by general duality theory 
that M' = IJpgf^'H-p, and we may equip M' with the inductive topology, 
that is the hnest locally convex topology such that all inclusions H-p C M' 
are continuous. We end up with the following chain of dense and continuous 
inclusions: 


M c Up+i cHpCHc U-p c 'H_(p+i) c M'. 


We also use tensor products and complexifications of these spaces. In the 
following we always identify / = [/i,/ 2 ] G "Hp^c,/i ,/2 G Up for p G Z with 
/ = /i + if 2 - The notation I'l^ is kept for the norm and (•, •) denotes the 
bilinear dual pairing. 

Remark 2.1. A/" is a perfect space, i.e. every bounded and closed set in M 
is compact. As a consequence strong and weak convergence coincide in both 
J\f and A/"', see page 73 in |GV64] and Section 1.6.3 and 1.6.4 in |GS68] . 

Example 2.2. Consider the white noise setting, where M = iS(M) is the 
space of Schwartz test functions, li = L^(M, dx) and M' = iS'(M) are the 
tempered distributions. iS(M) is dense in L^(M, dx) and can be represented 
as the projective limit of certain Hilbert spaces Up, p > 0, with norms 
denoted by I'l^, see e.g. |Kuo96] . Thus the white noise setting is an example 
for the nuclear triple described above. 

2.2 The Mittag-Leffler measure 

As Mittag-Leffler measures 0 < /9 < 1, we denote the probability mea¬ 
sures on M' whose characteristic functions are given via Mittag-Leffler func¬ 
tions. The Mittag-Leffler function was introduced by Gosta Mittag-Leffler in 
|ML05] and we also consider a generalization hrst appeared in |Wim05] . 

Definition 2.3. For 0 < /? < oo the Mittag-Leffler function is an entire 
function defined by its power series 
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Here T denotes the well-known Gamma function which is an extension of the 
factorial to complex numbers such that r(n-|-l) = n! for n G N. Furthermore 
we define for 7 G C the following entire function of Mittag-LefHer type, see 
also |EMOT55] . 

The Mittag-Leffler function is an entire function and thus absolutely con¬ 
vergent on compact sets. Hence we may interchange sum and derivative and 
calculate the derivative of E^: 


Lemma 2.4. For the derivative of the Mittag-LefHer function E^, > 0, it 

holds 



E/3,/3(^) 


^ G C. 


The mapping {fGM|f>0}9fi—)■ E^(—f) G M is completely monotonic 
for f3 G (0,1], see |Pol481 IFelTl] , This is sufficient to show that 


A/” 3 ^ I —V E^ 


;(e,e) e 


is a characteristic function on A/", see |Sch92) . Using the theorem of Bochner 
and Minlos, see e.g. |BK95] . the following probability measures on A/"', 
equipped with its cylindrical a-algebra, can be defined: 

Definition 2.5. For jS G (0,1] the Mittag-Leffler measure is defined to be 
the unique probability measure p /3 on A/"' such that for all G A/" 


J^exp{i{u,0)dix/3{u) = E^ ' 

The corresponding spaces of complex-valued functions are denoted by 
:= ^i 3 ]C) for p > 1 with corresponding norms By 

E^^(/) := /(cn) dfig^co) we define the expectation of / G 

In [Sch92] all moments of are calculated: 

Lemma 2.6. Let G A/" and n G N. Then 

[ and [ (a., 

JM' JN' ^ \P^ “T 

In particular = f(^ 1^1^- 

Remark 2.7. Using Lemma [2.61 it is possible to define (-jp) for rj & l-i as 
the L^(/r/ 3 )-limit of where is a sequence in M converging to rj 

in K. 
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Lemma 2.8. For ^,r] G Af it holds that 

Proof. By dehnition we have for each Oi, 02 € M 



exp {i{oj, ai^ + a 2 r])) = 




Applying the operator — to both sides and evaluate at ai = 02 = 0 

we get the result. Interchanging of sum and derivative is possible since the 
Mittag-Leffler function is entire. Further for all Oi, 02 G (—1,1) it holds 


d 

da I 


i{u, rj) exp (i(a;, ai^ + 0277 )) 


< |(a;,77)|el“2l<‘"’’?> |(a;,01 


g G L^{fii 3 ) by Lemma 4.1 in [GJRdSlH] . Hence d/dai interchanges with 
the integral. A similar argument yields that also the integral and 8/802 
interchange. □ 

Remark 2.9. It is shown in |GJRdS15] that / G Lf{gip) does not admit a 
chaos expansion if /? 7 ^ 1. This means that there is no system of polynomials 
(4(0) nSN’ eeAA, on M' fulhlling the following properties simultaneously: 

(i) IniO = Pn,0(‘4)) some monic polynomial Pn,^ on M or C of degree 

77. G N. 

(li) 4(0 T ImiO in 4(/i/3) for n ^ m. 

(m) 4(0 -L 4 ( 77 ) in 4(/i^) if ^ T 77 in 4. 

This is in contrast to the case /3 = 1 where the Mittag-Leffler measure is 
a Gaussian measure. It is a well-known fact in Gaussian analysis that such 
a system of orthogonal polynomials, the Wick-ordered polynomials, always 
exists, see e.g. |Kuo96[[HKF^ . 


2.3 Distributions and Donsker’s delta 

Instead of using a system of orthogonal polynomials, it was proposed in 
|GJRdS15] to use Appell systems, compare |KSWY98] . These are biorthog- 
onal systems allowing to construct a test function and a distribution space. 
As shown in |G,TRdSl^ the measures /r^, 0 < /9 < 1, satisfy the following: 
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(Al) The measure /i^ has an analytic Laplace transform in a neighborhood 
of zero, i.e. the mapping 

Afc3 := exp((a;, 9)) Q(6', 9)^ e C 

is holomorphic in a neighborhood U C Ac of zero. 

(A2) For any nonempty open subset U C Af^ it holds that fijsiU) > 0. 

We introduce the space of smooth polynomials on A/', denoted by V{Af') 
and consisting of finite linear combinations of functions of the form (-, 0 ”, 
where ^ G Ac and n G N. Every smooth polynomial (p has a representation 

N 

,.(0.) = 5^ (a.®”, 

where iV G N and is a finite sum of elements of the form ^ G Ac- 
We equip V{Af') with the natural topology such that the mapping 

CXD 

^ ,nen} 

n=Q 

becomes a topological isomorphism from V{Af') to the topological direct sum 
of tensor powers A"®", i.e. 


V(.M') ^ 0 V|“ 

n=0 

(note that 7 ^ 0 only for finitely many n E N). Then we introduce the 
space V'^^{Af') as the dual space of V{Af') with respect to i.e. 

V{Af') C C V'^^iAf') 

and the dual pairing ((•, •))^^ between V'^^{Af') and V{Af') is a bilinear ex¬ 
tension of the scalar product on by 

((/w))m^ = (/w)l2(^^), ip G P{Af')jE 

Note that (Al) ensures that V{Af') C L‘^{pp) is dense, see |Sko) . Then it is 
possible to construct an Appell System generated by the 

measure pp. We will only give the main definitions, for further details and 
proofs we refer to |G,lRdSL^ IKSWY9^ . 
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For each n G N and z G A/"^ we construct Pn^{z) G such that the 

Appell polynomials 


p/^/3 = |(P,^^(.),(pW) 
give a representation for ip G by 


^(n) e n G n| , 


AT 




( 1 ) 


71=0 


for suitable G N and G Every distribution <F G 'P^^(A/'') is of 

the form 

OO 

® = 5^ Qn" (2) 

71=0 

where is from the 



Furthermore the dual pairing of a distribution and a test function is given 
by 

{{Q'^rf ($(")) , = <5™,„n!(<F("),(p(")), n,m G No, (3) 

for G [Mc^^ G A/’^™'. With the help of the Appell system 

a test function and a distribution space can now be constructed, see 
[KSWY98], which results in the following chain of spaces 

C c L\pp) C c {M)-l p, g G N. 

Here, denotes the completion of V^M') with respect to || • 

given by 

N 

Ml.,.,, ■= P.9 e N„,v> e rm. 

71=0 

There are p',q' > 0 such that {'Hp)^ is topologically embedded in 

for all p > p',q > q', see |KK99] . Moreover, the set of /r^-exponentials 

|e^^(6 *;-) := \ d G Fp,g| is total in Here the neighborhood 









Up^q C Me of zero is defined by Up^q := G Me \ \0\p < 2 Since the 

exponentials have the series expansion ep^{9,-) = 
it holds that for each 6 G Up^ 2 q the norm is given by 

CXD 

n=0 

By {'H-p)z\^pp we denote the set of all <h G 'Pp^{M') for which 
is hnite, where 

OO 

ll'S’ll-™ := E^"’” P’1<=- No- 

n=0 

It holds that {l-L-p)z\^pp is the dual of {1-LpYqp^ and the test function space 
is defined as the projective limit of {'HpYq^^. This is a nuclear space 
which is continuously embedded in Moreover it turns out that the 

test function space is the same for all measures /i^ satisfying (Al) 

and (A2), thus we will just use the notation (A/")^. The space of distributions 
{M)~^ is the inductive limit of {l-L-p)z\^pp and is the dual of {MY with respect 
to M{^jY) and the dual pairing between a distribution $ G {M)~^ as in ([2]) 
with a test function tp G {MY as in ([T]) is given by 

OO 

n=0 

We shall use the same notation for the dual pairing between and 

_ 

As in |GJRdSl^ we introduce the -transform of d) G {M)^^ on a 
suitable neighborhood W C Me of zero by 

= ((<*.e„,(9, -)))„„ = 9 e M. 

Furthermore we define the T^^-transform of $ G {M)~^ at 0 G W by 

= ((*f’>exp(i(-,0)))^^. 

For <F G {M)~^ as in ([2]) we have 

OO 

= (5) 

n=0 
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The relation between and 5^^-transform is given by, see |GJRdS15] 

= eeu. (6) 

The space can be characterized via the -transform using spaces 

of holomorphic functions on A/c- By Holo(A/'c) denote the space of all 
holomorphic functions at zero. Let F and G be holomorphic on a neighbor¬ 
hood V,U <Z Me of zero, respectively. We identify F and G if there is a 
neighborhood W C V and W <Z U such that F{6) = G{9) for all 9 ^ W. 
Holo(A/c) is the union of the spaces 


F e Holo(A4:) 




7),/,oo 


{F) 


sup |T’(6')| < exo 
\e\p<2-i 


p, / e N 


and carries the inductive limit topology. 

The following theorem is proven in |KSWY98] : 

Theorem 2.10. The S^j^^-transform is a topological isomorphism from {M)~^ 

to Holo(A/'c). Moreover, if F e Holo(A/'c), F{9) = ^,11 

9 G A/c with \9\^ < 2“^ and if p' > p with < oo and g G N such that 

p ;= 2 ‘^’‘-^e^\\ipi< 1 then <l> = ^ 

As a corollary from the characterization theorem |GJRdS15] proves a 
result which describes the integrable mappings (in a weak sense) with values 

in 

Theorem 2.11. Let be a measure space and G (A/")”^^ for all 

t &T. Let U C Me be an appropriate neighborhood of zero and G < oo such 
that: 


(i) S^p^.{9): T ^ C is measurable for all 9 eU. 

(a) Jrp |S'^^d>t(6*) I dz/(f) < G for all 9 eU. 

Then there exists T G {M)~^ such that for all 9 eU 

Sw'I'W = «>.(«) dK*). 

ITe denote T by and call it the weak integral o/$. 

Another consequence of Theorem 12.101 characterizes the convergent se¬ 
quences in {M)~^: 
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Theorem 2.12. Let ($n)neN be a sequence in Then ($n)neN con¬ 
verges strongly in if and only if there exist p, g G N with the following 

two properties: 

(i) (S™-I>„(»))„a, is a Cauchy sequence for all 9 G Up^q. 

(a) Sfj,pC^n is holomorphic on Up^q and there is a constant C > 0 such that 

\S„iJ0)]<C 

for all 6 G Up^q and for all n G N. 

Proof. First assume that {i) and {ii) hold. From ([S]) and by Theorem 12.lUI it 
follows that there exist p', g' G N such that 

lltnil-p',-,',™ < np,,,oo(S™4„)(l - < C(1 - (7) 

Because of (i) together with ([Tj) and since the /X/^-exponentials are a total 
set in {PLpiYqi we conclude that the sequence (((<F„, is a Cauchy se¬ 

quence for all If G {Hpi)\i by approximating ip with /r^-exponentials. Indeed, 
for e > 0 choose e from the linear span of the /i^-exponentials such that 
\\p ~ cWp'^q'^pp < 1/(4C)(1 — pY^‘^e and choose m,n large enough such that 
i((<Fm - <hn,e))| < e/2. Then 

|{(*hn — l((*^n *h,T2,e))|-l-||*F^-|- || ||g? ® ^ 

Thus the mapping $: {UpYp p^ C, $((p) := lim„^oo((*hn, (^)), p G {Hpi)\, 
is well-defined, linear and continuous since 

|$(V3)| < \iuli■n^\\^n\\-p',-q',^,p\Mp'P,vp < p)-^/Y\t\\p’Pwp- 

n—^oo 

Hence <F G {l-i-pi)z\i This shows that (<hn)nGN converges weakly to $ in 
since {MY is reflexive. Finally we use that in the dual space of a 
nuclear space strong and weak convergence coincide, see Remark 12. 11 Hence 
converges strongly to <F in {M)p^. 

Conversely let (<Fn)„eN converge strongly to some $ G {M)~^. Strong con¬ 
vergence implies that there exist p, g G N such that <F in {l-L-p)z\^pp as 

n ^ oo. Thus (z) is obviously fulfilled for all 6 G Up^q. The strong convergence 
also implies that there exists K < oo such that sup^gj^ < K. 

Thus we have for all 6 G Up^ 2 q and for all n G N that 

KSp.tJWl < \\‘inU.-,„\\e„{ep)\\„„ < A'(l - 2-’)-*''^ < «D, 
see (jlj). This shows {ii). □ 
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Remark 2.13. Because of dH]) is holomorphic if and only is 

holomorphic. Thus the characterization theorems 12.101 12.111 and 12.121 also 
hold if the -transform is replaced by the T^^-transform. 

The characterization theorems have already successfully been applied 
in Mittag-Leffler analysis. A generalization of Donsker’s delta from Gaus¬ 
sian analysis could be constructed in via Theorem 12.111 and its T^^- 

transform is calculated, see |GJRdSl^ . In fact, for 6 TV and a G M 
Donsker’s delta is defined as a weak integral in the sense of Theorem 12.1 II bv 

= ^^exp(ia:((-,r/) - a)) dx G . ( 8 ) 

We prove in the following theorem that the weak integral in ([ 8 ]) can be 
approximate by a sequence of Bochner integrals from Thus Donsker’s 

delta can be represented as a limit of square integrable functions. 

Theorem 2.14. For rj G FL and a G M if holds that 
<^a((-,h)) = lim — / m 

n^oa ITI J 

Proof. Set := (27r)“^ n G N. First note that G 

L‘^{^jf) as a Bochner integral since 

/ n pn 

II exp (ia;((-, rj) - a)) ||l 2 (^^) dx= Idx = 2n < oo. 

■n J —n 

We have for the T^^-transform of that 

1 /I 1 \ 

see Proposition 5.6 in |GJRdSl^ . Now it holds that the integrand 
(a;)e"““E/ 3 (^-^a :2 ( 77 , 77 ) - ^{0,0) - x{e,r])^ 
converges pointwisely for each x G M to 

e-““E^(^_la;2(77, 77) _ 0 ) _ j 


as 77 —)■ cxD and it is bounded by 


E/3 (-|a;2(77,77) - 1 ( 6 , 9 ) - x(0,77) j 


It was 


shown in Proposition 5.2 in |GJRdS 15 ] that there is a neighborhood U C Me 


of zero and a constant G > 0 such that 

1 o, , 1 


E /3 (^--x^{ti, v) - 2 V) 


dx < G, 9 eU. 
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Applying dominated convergence we see that (9) converges as n —)■ 

oo to 

J M 

Moreover for all 6 ^ G W it holds 


(9)1 < 





di-<C. 


Now the assertion follows by applying Theorem 12.121 


□ 


3 Grey Noise Analysis 

3.1 Basic Definitions 

The main ideas of grey noise analysis go back to Schneider in |Sch90] . He 
constructed grey Brownian motion on a concrete probability space. Further 
details were given amongst others by Mura and Mainardi in [MMOQ] and also 
in |Kuo961 IKS93] . We now want to point out that grey noise analysis is a 
special case of Mittag-Leffler analysis, where the spaces W and "H are chosen 
in a suitable way as follows. 

Consider the space iS(M) of Schwartz test functions equipped with the 
following scalar product 

{^,v)a = C{a) [ \x\^~'^^{x)fi{x)dx, ^,r]eS(R). 

Jm 

Here, 0 < a < 2 and C{a) = r(Q; + 1) sin(^). The notation fj stands for the 
Fourier transform of r/ G iS(M), which is dehned by 

fj{x) = {Tr]){x) = / r]{t)e^*^dt, a: G M. 

\/27r Jr 

By II • II Q, we denote the norm coming from (•, Ja- Define the Hilbert space "Hq, 
to be the abstract completion of iS(M) with respect to (•, Jq,. |Sch92] gives 
an orthonormal system in iS(M) with respect to (•, •)« and |MM09] 

constructed an operator A" on "Hq, such that = ( 2 n + 2 + 1 — a)/i“. 

This allows to dehne the spaces 

:= {<^ G I IKATV^II, <cx)}. 

Remark 3.1. Dehne Sa := pr limp^^^ 'Ha,p and S'^ its dual space. The above 
considerations show that C Tia C 5^ is a nuclear triple. Thus we may 
introduce a Mittag-Leffler measure on S'^. Nevertheless we propose a 
slightly different choice of the nuclear triple J\f CH G M' for two reasons. 
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(z) The space l-ia is defined as an abstract completion. Therefore its el¬ 
ements are Canchy seqnences. For the definition of ggBm later on it 
is necessary to have l[o,t) G 'Hq,. The question arises in which sense 
the indicator functions l[a,b), a, & G M, are in Tia- A way out would 
be to identify the indicator functions with a sequence of Schwartz test 
functions (^n)neN) where (^n)nGN approximates the indicator function 
with respect to || • ||q,. This is still open. Further the later analysis 
becomes more complicated. 

(ii) The nature of the test function space S^, given as the projective limit 
of the spaces T-La,p is not obvious. It is not clear if Sa coincides with 
S{R). 

Therefore we find it useful to develop in the following an alternative approach 
to a grey noise analysis. 

We consider the usual nuclear triple from white noise analysis, see Ex¬ 
ample 12.21 The corresponding distribution space will be denoted by 

Define the operator on iS(M) by 

^ H = \, 

[Knlth, He 

with normalisation constant 

Remark 3.2. (z) and denote the right- and left-sided fractional 

derivative and fractional integral of order a, respectively. More details 
are given in Appendix below. Note that it is not important for 
the definition of , whether the fractional derivative is in the sense 
of Riemann-Liouville, Caputo or Marchaud, since all these fractional 
derivatives coincide on iS(M), see Theorem IA.161 below. A necessary 
condition for / to be in the domain of the Caputo derivative is that 
/ is differentiable. Hence the Caputo derivative of indicator functions 
is not defined. Therefore the operator denotes in the following 
the Riemann-Liouville fractional derivative or the Marchaud fractional 
derivative. 

(zz) Although defined on iS(M), the domain of is larger. In particular, 
can be applied to indicator functions l[a,b), —oo < a < b < oo, 
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and it holds for f G M that 


Im)) {t) = Yli + H-1/2) ^ '^0 ’ 

see Example I A. 71 below. Here and below, we use for a G M the con¬ 
vention 



f > 0, 
t < 0 , 


and t“ 


0, t > 0, 
(-f)“, t < 0 . 


Theorem IA.4I and Corollary IA.9I below show that 


e L^(M, dx). 


The following result for the Fourier transform of a fractional integral or 
a fractional derivative of 99 G iS(M) is known f |SKM93] ): 


= ip{x){Tix) " (9) 

TD'^(p{x)=(p{x){Tix)'^, (10) 

where 0 < a < 1, (=Fix)" = |a;|“ exp(=F^ sign(a;)) and cp G iS(M). With the 
help of this result we give a relation between the a-scalar product and the 
usual L^-scalar product, compare |LV14] . 

Lemma 3.3. The a-scalar product and the L^(M, da;)-scalar product are 
related by the operator . In fact for all ^,ri E iS(M) it holds 

Proof. From |Mis08] we know that the constant can be calculated to 
be Kff = ^y2Hsm{TrH)T{2H) and thus A '^/2 ~ C{a). Now let cp G iS(M). 

In the case 1 < a < 2 it holds that Mf/'^ip G L^(M, dx), see Theorem I A. 41 
below, since (p G iS(M) C L^(M, dx) for all p > 1. In the case 0 < a < 1 
note that dehnes a square integrable function. Thus (ITU]) 

and Plancherel’s theorem imply that G L^(M, dx). Hence, using 

([H]) and (ITU]) , we directly see that 

/ L2(R^dx) 

Now the assertion follows by Plancherel’s theorem. □ 
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Corollary 3.4. For every a, 6, c, d G M it holds 

(l[a,6), 1 m))„ = . 

Proof. From Remark 13.21 we already know that M^^‘^l[a,b) ^ F^(M, dx). For 
1 < a < 2 it is calculated in [SKM93] that for all ^ ^ (M, dx) it holds 

x G M. 

The same holds for 0 < a < 1, see Lemma fA. 101 below. The assertion follows 
now by the same arguments as in the proof of Lemma 13.31 □ 

Corollary 3.5. For all s,t >0 and for all 0 < a < 2 it holds that 

(aU-'U,,,,, = \{t° + S“ - l« - o|“). 

Proof. Use the previous lemma and Example 3.1 in [MMOO) . □ 

As in Dehnition 12.51 we dehne the corresponding Mittag-Leffler measures 
/i/ 3 , 0 < /5 < 1, on iS'(M) by 

f e‘<“’'">d,,^M = EA-ife4,)'), ^e5(K). 

JS'(R) V ^ / 

The measures /i /3 are referred to as grey noise reference measure. For /3 = 1 
this is the Gaussian white noise measure. Lemma 12.61 gives all moments of 
/i /3 and Remark 12.71 allows the extension of the dual pairing (•, •) to iS'(M) x 
L^(M, dx). Note that G L^(M, dx) for all f > 0 by Corollary IA.9I 

and Theorem IA.41 below. Thus the following dehnition makes sense: 

Definition 3.6. For 0 < a < 2 we dehne the process 

5'(M) 9 0;^ := {u, t > 0, 

and call this process generalized grey Brownian motion. In the case a = (3 
we write Bf instead of Bf is called grey Brownian motion. 

Remark 3.7. In the approach of |MM09] the grey noise measure is dehned 
via the characteristic function E/ 3 (—(•, •)q,) and denoted by /ia,/ 3 . This means 
that hrst the parameters 0 < a < 2 and 0 < /9 < 1 are hxed and then gener¬ 
alized grey Brownian motion is constructed in L‘^{na,/ 3 )- The measure 
/i /3 as dehned above is named grey noise reference measure since for hxed 
0 < /5 < 1 all generalized grey Brownian motions Bf’^ for 0 < a < 2 can be 
constructed in the single space 
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Proposition 3.8. Let 0 < a < 2 and 0 < /? < 1. Then for all p G N there 
exists K < oo such that 






<K\t-sr, 


t,s > 0 . 


Proof. Without loss of generality let s < t. First note that ||l[s,i)||Q, = 
||l[o,t-s)||a- Corollary 13.51 shows that ||l[s,i)|la = (f ~ Using the mo¬ 
ments of /i/ 3 , see Lemma [2.61 we find 


E 


'M/3 


B 






2p 


(2p)! 

2pr(/dp +1) 
(2p)! 

2pr(/5p +1) 


111 11 
(t - 5)“^* 


□ 


The last proposition ensures that each generalized grey Brownian motion 
has a continuous version. Indeed, choose p G N such that ap > 1 then the pre¬ 
vious proposition provides the estimate — Bf'^Y'P) < K\t — 

with q = ap — 1 > 0. This estimate is sufficient to apply Kolmogorov’s 
continuity theorem. 

Proposition 3.9. The processes [Bf’^] ,0<a<2, 0</5<l, have the 

V / i>0 

following properties: 

(i) has covariance 


E 




2T{I3 + 1 ) 


- |f - sD = 




t,s > 0 , 


and E,({Bf’^Y) = for f > 0. 

(ii) The density of the finite dimensional distributions of {Bf'^)t>o with 
respect to the Lebesgue measure are given by 


UA^) 


(2^)-n/2 p ^ I' 1 1 


where x G M” and ^ 0,0 denotes the matrix {pfaA'^i^'^j))i j=i n ^ 
For the definition of the M-function see Appendix [Bl below. 

(in) {BA)t>o has stationary increments. 
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Proof. The covariance can be calculated using Lemma 13.31 and Lemma 12.81 
and Example 3.1 in |MM09] . Obviously we can deduce that E — 


1 


r(/3+i) 


t". The second assertion follows by the same arguments as in |MPn8] . 


Proposition 1. Proposition 3.2 in |MMn9] and Lemma [3.31 show that 


E =E(^i 


iOB' 


. Q .,0 


This shows the third assertion. □ 

Remark 3.10. Consider the scaled ggBm := Then a compar¬ 

ison of Proposition 13.91 with Proposition 2 in |MPn8) shows that the scaled 
version of Definition 13.61 coincides with the dehnitions of |MPn81 IMMn9] . 
Consequently B^’^ is a grey Brownian motion as in |Sch92] . 

Remark 3.11. Generalized grey Brownian motion is a generalization of well 
known stochastic processes. Choosing (3 = 1 then Proposition 13.91 shows 
that E^j = l/2((t" -|- s" — |t — s|“) and the hnite dimensional 

distributions of Bf’^ are given by 


fa,i{,x) 


1 

— , exp 

^27rdet7a,i 



Thus Bf'^ is a fractional Brownian motion with Hurst parameter H = a/2. 
Consequently the process B]'^ is a Brownian motion. 

Remark 3.12. A random vector X = (Xi,..., X„)^, n G N, is said to have 
an elliptical distribution with parameter m G symmetric and positive 
dehnite S G and some characteristic function 0: M —)■ M (for short 

X ~ EC„(m, S, (/>)) if its characteristic function is of the form 


E(e‘®^^) = e G 


see e.g. equation (2.11) in |FKN90] . The proof of Proposition 13.91 shows that 
X = ... ,Bff) has an elliptical distribution and X ~ EC„(0,r(/5 -|- 

1 )7 q,^^,E/3 (—)), see also |dSE15] . 


4 The time-fractional heat equation 

The well-known heat equation 

d 

=2 

M(0,a;) = uo{x), 


( 11 ) 
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for t > 0, a; G M, and initial value Uq G iS(M) can be generalized to a 
fractional differential equation in two ways. Replacing the time derivative by 
the Caputo derivative of order a denoted by ‘"£* 0 +) 0 < a < 1, we obtain the 
equation 

1 8 “^ 

{^DQ^u){t,x) =-—u{t,x), xER,t>0 . (12) 

Starting from the heat equation in integral from and replacing the integral 
by a Riemann-Liouville fractional integral gives the equation 

1 /■* d'^ 

u{t,x)=uo{x) + -— {t - s)°'~^—u{s,x)ds. (13) 

r(a) Jq dx^ 

flT^ has been solved for example in |Mai95] and a solution to ffT3|) is given 
in |SW89] . The Green’s function G has the series expansion 


G(t, X, a) = 


\/4^fro 


E 


(-1)^ T{l/2-k) 


k\ r(l — ak 




E 

/c=0 


(- 1 )^ r(-l/ 2 -A;) 
k\ r(l — ak — a) 


r 

a/2) 

9 \ fc 


(14) 


The same solution has also been given in |Koc90] . Note that in general 
fractional differential equations with Caputo derivative are equivalent to the 
fractional integral equations using the Riemann-Liouville fractional integral 
if the solution is absolutely continuous, see Theorem IA.121 below. 

For the heat equation (ITT]) the solution is given by the Feynman-Kac 
formula, i.e. u{t,x) = E(mo(3^ + Bt)), for twice differentiable initial value 
Mo with compact support and a Brownian motion B. |Sch92] showed that 
u{t,x) = K{uo{x+B^)) solves the fractional heat equation flT^ for uq G iS(M) 
and grey Brownian motion '■= {Bnt>o. We prove in the following 

that, replacing uq by 6^, a Green’s function to flT^ is obtained. For the proof 
we need the following lemma: 

Lemma 4.1. For all A G C and for all (3 G (0,1) the Mittag-Leffler function 
Ejs satishes the following equation: 


Proof. First a simple calculation shows that 


a/l3-l 





ds. 


ds. 


(15) 
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Indeed, insert the series expansion of in the right hand side. Due to the 
absolute convergence of the series we may interchange sum and integral. For 
the general case a ^ (3 note that 


Applying flT^ yields 


E/3 



1 - 


1 \ 2 ^_ 

2 r(/3) 





ds. 


(16) 


Now the coordinate transform s = yields the desired result. □ 

Remark 4.2. Equation flTB]) shows that 




(17) 


Theorem 4.3. Let 0 < /9 < 1 and 0 < a < 2. For x, ?/ G M and f > 0 
define K{t, x, y) := {^y{x + . Then K is a Green’s function to the 

equation 

u{t,x) = Uo{x) + ^ f > 0, X e M, 

with initial value uq G L^(M, dx) Pi L^(M, dx) satisfying 


^mo)(A)A|^'''^ dA < cx) and I |(J^ ^Mo)(A)A^h'^^ dA < oo, 


I l+£ 


for some £ > 0. 

Proof. Define u(t,x) = J^uo{y)K{t,x,y) dy for x G M and f > 0. Note that 


E 






see Proposition 5.6 in |GJRdS15] . Equation (l3^ in the appendix below 
implies that 

uo{y)^up {^yi^ + Mg{r)r-^/^drdy. 
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Since Mi 3 {r)r dr = vr/r(l — a 




dy < 


Hence u is well defined and it holds 


/2), see e.g. |MMP10] . we obtain 
,_ —^ - / Inofy)! dy < oo 


u{t, x) = 







uo{r)e 



dA 


Denote the integrand by f{x, A), x, A G M. Then f{x, •) G L^{R) for all x G M 
and /(•, A) is differentiable for all A G M with 


A 

dx 


f{x,X) 


|A| 


E. --XH' 


2.a 


B-\o) (A) 


:= ^(A). 


We show that g G L^(M, dx). With equation 0351) below we estimate g as 
follows: 


5 '(A)dA< / Mjs{r) / ^mo)(A)A exp —-A^f"r dAdr. 


Next we use Holder’s inequality. Choose p = 1 + 1/s and q = p/(p—l) = l+e. 
Then l/p + l/g = l and 


|(J^ ^Mo)(A)A|exp ( —-A^f“r ) dA 


< 


1 

(J'”^Mo)(A)Ap’^^ dA^ exp (^~pX‘^t°‘r^ dA^ 


(J^-'no)(A)A|^+" dA 


TTi 


.-l/(2p) 


P+ 


Due to the condition on uo together with Lemma A.2 in |G,TRdS15] it follows 
that g G L^(M, dx) and we may interchange derivative and integral and obtain 
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With similar arguments we can show that also the second derivative inter¬ 
changes with the integral. Thus 

The assertion follows by using Lemma 14.11 that □ 

Remark 4.4. The series expansion of u for the case a = f3 is calculated in 
Appendix iBl below, as: 


u{t, x) 


1 A (- 1 )^ r(i/2-fc) 

\/27rf" “ k\ T{l-l/2a-ak) 


+ 


X 


OO 


(-1)^ r(-i/2-fc) 


r(i 


a 


ak) 



t > 0, X G M. 


Comparing this to (ITT|) we see that our approach is equivalent to |Mai95| . 
Remark 4.5. Dehne for f > 0 and x G M the function u{t,x) = E^^(5(x -|- 
Then it can be shown that u solves the equation 


u{t,x) = So{x) + - 


r(/?) dx^ 


m(s^/“,x) ds 


in the distribution space iS'(M). Note that 


u{t, x) = 


1 dA 


27r 


2 +a 


Inserting the resnlt from Lemma 14.11 we achieve: 




= — / e'^" dA 

StT /to 


y 2 


r(/9) 


E, 


--\h^ 1 ds 1 dA 


1 (^a//? _ 5)/3-1 Q2 / 1 


— ho(3^) + X 


r(/ 3 ) dx'^ \2 'k 


dA ds 


r(/?) 9x2 


x) ds. 


The calculation and the involving divergent integrals can be treated rigor¬ 
ously by applying them to a test function ^ G iS(M). 
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Remark 4.6. For the special case of fractional Brownian motion we set 
/3 = 1. Then we obtain that 


K{t,x,y) = {6y{x + , t > 0, x,y eR, 


is a solution to 


u{t,x) = Uo{x) + - 

= +1 


I 


—ds 


I 


as" M(s,x)ds, 


or equivalently 



(19) 


In the case of the Gaussian measure /ii the T^^-transform of Donsker’s delta 
can be calculated explicitly: 



It can be easily verihed that K is indeed a solution to flT^ . 

5 Grey noise process 

We have seen in the previous section that generalized grey Brownian motion 
plays a central role in the context of the time-fractional heat equation. 
This motivates the further analysis in the following section. We hrst calculate 
the 5^^-transform of Then we prove that the local times of generalized 

grey Brownian motion exist. 

In the following we calculate the -transform of a generalized grey 

Brownian motion = (•, G for (p = (pi + iip 2 G W/? C 

iSc(M). Here, = {(p G Me \ |(<p,(p)| < Sy}, where > 0 is chosen such 
that (z) > 0 for all \z\ < By, z E C. Then the 5^^-transform of is 
well-defined and it holds that 


(V) 


E/3 


1 
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For u G iS'(M) and s G [—1,1] denote f{uj,s) := exp((ci;, (^ + 

Then / is integrable with respect to u and differentiable with respect to s 
and 

Moreover we have for all s G [—1,1] and for all u G iS'(M) the estimate 


d 



< 


{u!,MJ l[0,p) 

< e<‘"’^i>exp (2 |(a;,M“/2l[o,p)|) =: g{uj). 

It is shown in Lemma 4.1 in |GJRdS15] that both mappings exp((-, tpi)) and 
exp ^2 (•, j are in G Hence we have that g G L^{g>/ 3 ) by 

the Holder ineqnality. Thns interchanging derivative and integral is possible. 
Since satisfies the assnmption (Al) on page E] we get 

(^) 

see Corollary 4.3 in |GJRdS15] . Using Lemma 12.41 and evalnating at s = 0 
we obtain 

/' / /\W 2 i \ V^)) 

S„B, = 


We may use the integration by parts formula (12^ in the appendix below 
in the case 1 < a < 2 since ip and l[o,t) are in every L^(M, dx), p > 1. In 
the case 0 < 1 < a we are allowed to use fl27|) in the appendix below, since 


(i-G/2, 


p G L^(M, dx) and G L^(M, dx) for p = 2/(1 + 2a) by 


D 

Corollary IA.91 below. Thus we get 


(v) = l[o,.)> 


a/2 


(3E^{l{p,p)) 


) (x) dx. 


/3E^(l((p,(p)) Jo 

Remark 5.1. In the white noise case a = /3 = 1 we obtain 

l[o,t)) = / da; 

Jo 

which is the F-transform of Brownian motion. 


( 20 ) 
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Theorem 5.2. Each generalized grey Brownian motion (0, cxd) 3 t G 

for 0<a<2, 0</9<l, is differentiable, i.e. 


lim 

h^O 


B 


a,/3 

t-\-h 


B 


a,g 


h 


exists and converges to some element in 


denoted by 


fulfilling 


{s„Nrfi M = (m:'v) (i) 


E/3,/3 

f3Eg{^{(p,(p)) ’ 


for every cp from a suitable neighborhood U C Me of zero. 


Proof. Consider for n G N and t > 0 


:= 


B 


q,/3 

t-\-hn 


B 


a,y 


hn 




for a sequence (/i„)neN such that hn —t 0. By (120]) 


= T- 


f*t-\-hn 


h 


n Jt 


] {x) dx 


E/3,/3 (!(</?,</?)) 

fiEp{l{p,p)) ’ 

for every p G lAg. Since is continuous, see |Ben 03], we have 




lim -— 

n^oo h. 


n Jt 


) (x) dx 


E/3,/3 {\{p,p)) 
l3Eg{^{p,p)) 


/2,f\ 

l3Eg{l{p,p)) ■ 


= Ut) 


Thus {S^^^n(,p))n£N E & Cauchy sequence. Since the Mittag-Leffler functions 
are holomorphic there are p, g G N such that 


Eg^g{^{p,p)) 

fiEg{l{p,p)) 


for all p G Up^q. 


Furthermore it follows from Theorem 2.3 in |Benn3) that there is p' G N such 
that for all X G M 


( 3 ^) <C\p\ 


Here, (|•|p)pGN denotes the family of norms which defines the topology of 
iS(M), see Example 12.21 Choose p* > ma.x{p,p'). By the mean value theorem 
hnd s G [t,t + h] such that for all p G Up*^q and for all n G N 


s„ (B,th - Bp) 

M) 

< 

E/3,/3 {\{P,P)) 

hn 

l3Eg{\{p,p)) 




< KC 


\p' 
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Finally, applying Theorem 12.121 we see that (*hn)„gN converges to the distri¬ 
bution in {S)~^ fulfilling 

(^) = lim (f e Up*^g. 

\ / n—^oo 

□ 


5.1 Local times of generalized grey Brownian motion 

The ggBm local time, denoted by La,i 3 {a,T), measures the time for 
spending up to time T > 0 in a G M. It is already shown in |dSE15| that 
generalized grey Brownian motion admits a square integrable local time by 
checking Berman’s criteria. In the paper |dSE15j the local time is introduced 
as the Radon-Nikodym derivative of the occupation measure. In detail: Let 
/ be a measurable set in the interval [0, T], T > 0, and /: / ^ M measurable. 
The occupation measure /x/ is defined as 

If p/ is absolutely continuous with respect to the Lebesgue measure dx we 
denote the corresponding density by L-^{■,!), i.e. 

^ 1b(/(s)) ds = fif{B) = L^(x, I) dx, B G B(M). 

is called local time of /. Informally the local time can be expressed 
via the Dirac delta distribution. In fact it holds that 

Lf{xJ) = J^6,{f{s))ds, ( 21 ) 

see e.g. |Kuo96| . page 273. Motivated by (ET]) we give a second approach for 
constructing the ggBm local time T). For T > 0 and a G M we define 

Wa.T) := / (5„(B?-'’)d(€(S),7‘ 

Jo 

Theorem 5.3. The local time La,i 3 {a,T) of generalized grey Brownian mo¬ 
tion with d<a<2, d<l3<l,d<T <oo and a G M exists as a weak 
integral in {B)~^ with T^^-transform 

(T„L^_f(a,T)) K) = / K)d( 

for all f from a suitable neighborhood U C iSc(M) of zero. 
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Proof. Let ^ G >Sc(M) with |'ClL 2 (iRdx) ^ ^ some M < oo. According to 
Proposition 5.6 in |G,TRdS15] with rj = G L^(M, dx) we have 


(?) 

exp(-ixa)E^ - ^{?.?) “ a^(?. V^njo,,)) j dx. 


1 


The same calculation as in the proof of Proposition 5.2 in |GJRdS15] shows 
that 


T,MBP(0 


dt 


< 




(*T poo 


0 ^0 


M^(r)r"^/^f"/2exp drdt, 


with z{t,^) = (^ 1 , 6 ) - ( 6 , 6 ) - > -M^ since | 6 L 2 (R,da,) < ^ 


and thus 


TMBpm 


dt 


< 


\/^2 — a 


yi-a/2 / / ijVfV ) dr. 


By Lemma A.2 in |GJRdS15] we get 


T,MBpm 


di < 

V^2 - a 


for some constant K < oo. The assertion follows by Theorem 12.111 
Remark 5.4. The expectation of La_^(0,T) can be calculated explicitly: 


□ 


E(L„,^(T,0))= / E{6o{Bf’^))dt 


■ dt = 


6^r(l-/3/2) T{l-l3/2)y/^2-a 


rjnl—a/2 


This result extends the existence result for ggBm local time in |dSE15] . 
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A Basics of fractional calculus 


In this chapter we introduce the basic notations of fractional calculus, for 
a detailed overview and proofs we refer to |SKM93] . First of all we give 
the dehnition of the Rieniann-Liouville integral and derivative. Besides that 
there are a couple of other dehnitions. Mostly common in physics is the 
fractional derivative in the sense of Caputo. We also consider the Marchaud 
fractional derivative. This derivative appears in the context of fractional 
Brownian motion in the white noise setting. In general, these three fractional 
derivatives are different. But we prove that they coincide on iS(M). 


A.l Riemann-Liouville fractional integral and deriva¬ 
tive 


We hrst give the dehnitions of the fractional integral and derivative in the 
sense of Riemann-Liouville for functions / which are dehned on an interval 
[a, h] for —oo < a < h < oo. 

Theorem A.l f |SKM93] . Theorem 3.5). The Riemann-Liouville fractional 
integrals of order 0 < a < 1 of f E {[a, b], dx), defined by 


(C/) : 

1 

^r(a) 

J a 

— t)°‘ ^ dt, X E [a, b], 

(22) 

(lU) (x) : 

1 

^ r(a) 

f fm- 

J X 

- x)°‘~^ dt, X E [a, 6], 

m) 


exist if 1 < p < 1/a. Moreover /"+/ G T^([a, 6 ], dx), where 1 < q < U 

1 < p < 1/a then lQ_^_f G L'^{[a,b],dx) forq = fl22ip are called 

left-sided and right-sided Riemann-Liouville fractional integral, respectively. 

Theorem A.2 f |SKM93] . Lemma 2.2). Let f: [a,b] —)■ M be absolutely 
continuous. Then the left-sided (right-sided) Riemann-Liouville fractional 
derivative of order 0 < a < 1 of f given by 


{D‘j) (n 
(Dt-f) (x) 


r(l — a) dx Ja 
-1 d r’^ 
r(l - a) dx X 


t) “ df, X E [a, b], 

(23) 

x)~°' df, X E [a, 6 ], 

l|23I) 
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exists. Moreover D^_f G LP{[a, b] dx) for 1 < p < 1/a and 


(^5."+/) M = ((^ +1 /'(*)(- - , 

(«?-/) (-) = ru^ if3r - ■ 

Remark A.3. For functions / defined on the positive real line, the defini¬ 
tions fl 22 |) and fl 2 ^ are valid for all a: G [ 0 , cx)). 

Next we consider functions / which are dehned on the real line. 

Theorem A.4 f |SKM93] . Theorem 5.3). For every f G L^(M, dx) with 1 < 
p < 1/a the left-sided (right-sided) Riemann-Liouville fractional integral of 
order 0 < a < 1 zs defined as follows: 

(Ilf) (a:) := r f{t)(x - tr-' dt, I e R, (24) 

r(a) 

{Iff) {x) := / f{t){t - xr-^ dt, xeR. m) 

r(a) 

Moreover /^: L^(M, dx) —>■ dx) is continuous for 1 < p < 1/a and 

Q=t^- 

^ 1—otp 

In order to dehne the fractional derivatives we need the following: 

Lemma A.5 f |Mis08] . Lemma 1.1.2). Let 1 < p < 1/a and assume that 
/ G I^{LP{R,dx)). Then there is a unique p G LP(R,dx) such that f{x) = 
{i±f) (x). This <p is called Riemann-Liouville fractional derivative of /, 
denoted by Dff and given by: 


{Dtf) (^) 
(Dff) (^) 


1 d r 

r(l - a) dx 
-1 d /■“ 

r(l - a) dx X 


f{t){x-t) "df, 
/(f)(f-x)“" dt, 


X G M, (25) 

X G M. ([251) 


Remark A.6. (z) By the dehnition of the Riemann-Liouville fractional 

derivative it holds automatically that 


T>±/±/ = f, for all / G L^iR, dx), 1 < p < 1/a. 
Moreover, if / G/^ {Lp(R, dx)), 1 < p < 1/a, then we have also 

IfD^f = /• 
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{a) Equation (5.16) in |SKM93] proves the following integration-by-parts 
formula: Let / G L^(M, dx) and g G L^(M, dx) with p,q > 1 and 
- = 1 -|- a. Then it holds for all 0 < a < 1: 

p ' q 

[ (/+/) {x)g{x)dx = [ f{x) {^g) (x) dx. (26) 

M. v/ M 


[iii) The next formula gives the integration-by-parts formula for the frac¬ 
tional derivative, see (5.17) in |SKM93] : Let 0 < a < 1, T*"/ G 
LP{R,dx), D‘^g G L^(M,dx) with l/p-|-l/g = l-t-a and / G L'’(M,dx), 
g G L^(M, dx) with 1/r = 1/p + a, 1/s = 1/q — a. Then, 

[ (T>“/) (x)^(x) dx = [ f{x) (Dig) {x)dx. (27) 

t/ M t/ M 

As an example we state the fractional integral and derivative of the indi¬ 
cator function l[a,b), see Lemma 1.1.3 in |Misn8] : 

Example A.7. The indicator function l(a,fe) for a, 5 G M lies in the domain 
of /^. Moreover 

(^±l(a,6)) (t) = ± (« - t)^) • 

The Riemann-Liouville fractional derivative of the indicator function l(a,b) is 
given by 

(^±l(a,f>)) {t) = {T{b - t)-^ ± (a - t)-^) . 

Lemma A.8. For 0 < a < 1 it holds that 71±l[a,6) G L^(M, dx). 

Proof. See also |ST94] . For x < a and for some G [a, 6] we have by the 
mean value theorem 


(6 — x)_,_“ — (a — x)_,_"| = {h — a)a{^ — x) “ ^ < (& — a)a{a — x) “ 


Thus 



/ a—1 

(a — x)“"“^ dx < cxo. 

■OO 


For X G [a — l,a-|-(5], 5 > 0 small enough such that a + 5 < h., there is a 
constant C < oo such that {h — x)~°‘ < C and we get 



(a — x)_,_“| dx < C(1-|-5)-t- f (a — x) "dx<cx). 

Ja-l 
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For X G [a + 5, 6], we only have to consider the first term and find that 


{b — x) " dx < cxo. 


' a+5 


For X > b the function is zero. This shows that D^l\a,b) £ T^(M, dx). Use 
similar arguments to get £ L^(M.,dx). □ 

Corollary A.9. A slight modification of the previous proof shows that 
D±'dla,b) £ LP{R,dx) whenever ap < 1. 

Lemma A.10. Let 0 < a < 1/2. Then for every a, 6 G M it holds 


D^l[a,b]{,x) = (=Fia;)“l[a,6](a;). 

Proof. First let 0 < a < 1. From |GS64] it is known that a Fourier transform 
pair is given in the distributional sense by 


t 


0—1 

± 


r 
i — 


■» 



{Tix) “. 


This means that for each p G iS(M) it holds that 


y In 


i.e. 

[ = j {TixY~^p{x)dx. 

Jr V Jtt Jm 

Set f{x) := — (6 — a;)l" + (a — a; G M. Then Corollary IA.9I shows that 
/ G L^(M, dx) n L‘^(K,dx) for 0 < a < 1/2. Furthermore it holds for each 
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if G iS(M): 


= / f{x){J^^p){x)dx 


= — {b — x)_°'{J^ip){x) dx + / {a — x)_'^{J^ip){x) dx 

tJ M. t/ M 

POO POO 

= — {x — b)~°‘{J'Lp){x) dx + / {x — a)~°‘{J'ip){x) dx 

Jb J a 

POO POO 

= — t~°‘{J'ip){b + t) dt + / + t) dt 

Jo Jo 

= [ + - {J^ip){b+ t)) dt 

Jr 
1 



r(l -g) 
r(l-g) 


ds dt = / ((e“- - (^) dt 

Jr 

[ (-ia;)“-i(e“" -e''’")v9(a:) dx 

Jr 


Since iJ"l[a^fe](a;) = see Example I A. 71 we conclude that 

1 


Dll[a,b]{x) = -=(-ix)" (e^ 




a—lf^ibx _ 


Since it is well known that 

we obtain 


^ibx _ ^iax 

\P1t: ix 


D^l[a,b]{x) = l[a,b]{x){-ix){-ix)'^ ^ = l[a,b)(a?)(“io;)". 
A similar calculation for g{x) = (6 — x)+“ — (a — a;)+" shows that 

U"l[a,6](a;) = {ix)°‘lia,b]{x). 


□ 


A.2 Fractional derivative of Caputo-type 

In contrast to the Riemann-Liouville fractional derivatives (123]), (ESH) and 
(I25D, (EHj) we denote the fractional derivative of Caputo-type by '"ZJ". 
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Theorem A.11 ( |KST06] . Theorem 2.1). A sufficient condition for the left- 
and right-sided Caputo fractional derivatives of f : [a,b] — )■ M of order 0 < 
a < 1 is that f is absolutely continuous. Then 

(^T>“+/) (a:) := (T>“+/) (x) - 2 : G [a, 6], (28) 

(%“_/) (x) := {Df_f) (x) - X e [a, b]. ([28l) 

In this case it holds that (125]) and (® coincide with the following expres¬ 

sions: 

i^Df^f) {x) = f'{t){x -t)-^dt, X e [a,b], (29) 

i^D^-f) (x) = ^ f'{t){t - x)"“ dt, a; e [a, b]. ([291) 

Theorem A.12 ((D ie IQ], Theorem 3.7 and 3.8). Let f be continuous on 

[a, b]. Then 

= /. 

If f is absolutely continuous, then 

(/“+^T>“+/) (x) = f{x)-f {a), xe [a, b]. 

For /: M —>■ M the Caputo fractional derivatives on the real line are 
defined as follows, provided they exist: 

(Wlf) ]_ . f - t)~° it, a:eK, (30) 

— 1 

{Wlf)(x) = — - - f'(t){t-x)-“dt, leR. BH) 

r(i-a)7^ 

A sufficient condition for the existence of^DI f is that f G iS(M), see Theorem 
[AtHI below. 


A.3 Marchaud fractional derivative 


Next we define the left- and right-sided Marchaud fractional derivatives on 
the real line, provided they exist, for / : M —)■ M by: 


^’Kf) (X) ^ ^ 

L 


— ^de. 

a: G M, 

(31) 

^D^f) (x) - „ “ , 

roo 

f{x) - f{x 

+ de. 

a: G M. 

CD) 
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We also introduce the truncated Marchaud fractional derivatives on the real 
line by 


("o;,/) ix) = de, ^eR. m) 

Then 

("C”/) = lim ('•‘D% J) , 

where the limit is dehned by the problem under consideration. For /: [a, b] —)• 
M the Marchaud fractional derivative on the interval [a, b] is dehned for x G 
[a, b] by 


rDu) (-) = + 


a 


_ r fix)-fix-0 

r(l - a) A 

« f’’ fix) - fix - 0 


de, (33) 

de. dMi) 


r(i - «) r(i - a) A e+' 

Theorem A.13 f |SKM93] . Theorem 13.1.). Let f G LP{[a,b],dx). Then 


Mj~\OL ja £ _ £ 

^a+-^a+J ~ J • 


In particular, the Marchaud fractional derivative and Riemann-Liouville frac¬ 
tional derivative on an interval coincide for f G (LP([a, 6], dx)). This is 
true, if e.g. f is absolutely continuous. 

Theorem A. 14 f |SKM93] . Theorem 6.1). If f E dx) for 1 < p < 1/a. 

Thcfi 

(MDiiy) (j,) = f(x). 

Here ^Dff is the L^-limit of the truncated Marchaud fractional derivative 

m and dHl). 

Example A. 15. The Marchaud fractional derivative of the indicator func¬ 
tion l[a,6) is given by: 

("DJ (f) = (t( 6 - ± (a - t)-p . 

Proof. 11 X > b, then clearly i^Dfl^afi)) ( x) = 0. Let now x < a. Then 


(“flUla.i)) (X) = 


a 


l[a,6)(2: + 0 


a 


r(i - a) 


r(l - a) it 

r“"'de = 






r(l-a) 


((6-x)-“ - (a-x)-“) . 
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For X G [a, b] it holds 


{^D^.l[a,b)) {X) = 


1 - l[a,fe)(a; + 0 


r(i - a) 


r(l - a) 7o 

OO -1 

A b-x)- 


r(l - a) 


Altogether we obtain 




A similar calculation shows 


'B:i|a,«) (i) = (-(6 - i):“ + (a - *):") ■ 


This hnishes the proof. 


A.4 Fractional derivatives of Schwartz test functions 

Theorem A.16. The Riemann-Liouville fractional derivative coincides with 
the Marchaud fractional derivative and with the Caputo derivative on S (M). 

Proof. The integral f(xT ds for / G iS(M) exists since 


/ \f{xTs)\s “ds= / \f{xTs)\s “ds+ / \f{xTs)\s "ds 
Jo Jo Jl 

P 1 i*00 

< ll/lloo / s“"ds + sup |s/(s)| / < CX). 

Jo SSR Jl 

The coordinate transformation s = x — t yields 

= -1)-<-it 

= r,.^ . f- r Ax-,)s-'--ds. 
r(l - a) dx Jq 

Interchange of derivative and integral is justihed since \f'{x =F s)| is dom¬ 
inated by the integrable function 

h{s) = ||/'|lool[ 0 ,l)(s)s““ -t- sup |s/'(s)| l[i,oo)(s)s"“"^ 
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This yields 


iKf) W = 


r(l - a) Jo 

By similar arguments we obtain 


(D“f) (x) = 


-1 


f'{x — s)s " ds. 


f'{x + s)s “ ds. 


r(l - a) Jo 

Comparing with flHOj) and 03011) it follows that for all / G iS(M) 

D^f = ^Dlf. 

Furthermore since s“" = a d.^: 


(^:/) (^) = 


a 


f'(x-s) / ^ " M^ds 


r(l-a)io ./s 

r r/'(x-s)r“-Msde 

r(l-a)io Jo 

a fix)- fix-^) 


r(l - a) Jo 


e 


Q + l 


■de 


And analogously 

(oy) (^) = 

and 03111) shows that 


a 


r(l - a) Jo 


fix) - /(x + 0 


de 


Dlf = ^Dlf, feSiR). 


□ 


B The H-function 


The H-function was discovered by Charles Fox in 1961 |Fox61] and is a gen¬ 
eralization of the G-function of Meijer. The definition is as follows: Let 
m,n,p,q & N, 0 < n < p and 1 < m < g. Let Ai,Bj G M be positive and 
Oj, bj G M or C arbitrary for 1 < i < p, 1 < j < g. Then 


Tjm n 



(op, Ap) \ 

KB,)) 



(®i) ^i)) • • • 1 i^pi -^p) ^ 
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where 



(n", r(% + sB,)') (n”., r(l - aj - sA,)) 


For further details concerning the contour £ and existence of If we refer to 
[MSHIO] , The series expansion of if for \z\ > 0 is given in |PBM90] 



under the condition that > 0 Bk{hj + 0 7^ Bj{hk + s) 

for I < j,k < m,j ^ k and /, s G Nq. 

Remark B.l. A change of variables in the definition of the if-function 
shows the following properties, see also e.g. IKSTOh] : 

(i) For all z G C \ {0} such that the ii-function is defined the following 
inversion formula holds; 



(u) For all (T G C and for all z G C such that the if-function is defined 
there holds the formula 



The M-Wright function Mp for 0 < /9 < 1 was introduced by Mainardi as 
an auxiliary function when finding the Green’s function to the time-fractional 
diffusion-wave equation, see e.g. |Mai96bl IMai96a] . Its series expansion is 
given by 



For all 2 ; G C the following holds |GJRdS15] : 



(35) 
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For more details we refer to |MMP10] and the references therein. 
Let us now compute the expectation of grey Donsker’s delta: 


E 




a + Bf)) = (r,./(a + Bf)) (0) 

exp(ia;(5f + a))] (0) dx 


1 




27r 


Ma(s)exp ( dsdx 


Mi 3 {s) [ exp (--x^t ^s + ixa\ dxds 


27r 




Mp{s)s ^/^expf-^s ds. 


We express the Mainardi function and the exponential by H-functions. First 
we have 


= Hl« ( s 

Thus we get 


(1-/3./3) 

( 0 . 1 ) 


and e ^ = HH ( 2 ; 


( 0 , 1 ) 


E^, idia + B^) = 


V2Trt^ 




(l-A/S) 

( 0 , 1 ) 

^10 f 
»' Utfl 


( 0 , 1 ) 


ds. 


Using the inversion formula for the H-function, see Remark IB. II we end up 
with 


E„pi(o + sf)) = 


1 


\/2'KtP 


( s 


(l-/3,/3) 

( 0 , 1 ) 

( 1 , 1 ) 


f2t^ 

rOl / 


Lemma 2 in |MS69] shows that 


cr—1 ijmn 
PQ 




ax 




{Cry Cr) 
{ds,Ds) 


ds. 


dx 


= a-’H^tPitT (ff/a 


{1-bq- Bqa, Bq), (c^, Cr) 

(1 -d-pC, -dp), (ds, Dg) 
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under the following conditions 


3?((j + bj/Bj + di/Di) >0, j = 1,... ,m, i = 1,..., k, 

3?((T + {uj - l)/Aj + (cj - 1)/Cj) < 0, j = 1,... ,n, i = 1,..., I, 

m q n p 

-^1 = X] ^ ^ + X] ^ > 0’ 

j=l j=m+l j=l j=n+l 

k s I r 


^2 = E 

i=i 


E Di + Y^Ci 

j=k+l j=l 


5^ Q > 0, 

j=l+l 


arga| < Ai7r/2, |arg/9| < 


In our case we have a = 1/2 and m = p = q' = l,n = 0 and a = 1 and 
k = s = 0, l = r = l and {3 = 2t^ /a? and Ai = 1 — /3 > 0 and A 2 = 1 > 0. 
Thus all conditions are fulhlled and we obtain 


(i(a + Bfi) 




tj02 

-"21 



(1/2,1),(1,1)A 

(/9/2,/9) )■ 


By the inversion formula from Remark IB. II it holds 


E;., (5a((-,h))) 


^ jj20 f 
^^2'K{T],r]) V2(h,h) 


(l-/5/2,/9) \ 
(1/2,1), (0,1);- 


The dehnition of the TT-function implies 


Em, (^a((-,h))) 


^ ^20 f 

^j2'K{ri,ri) \2{rj,vi) 


(l-/3/2,/3) A 
(0,1), (1/2,1);- 


Applying fl3T)) we hnally have 


Em, (^a((-,r;))) 


1 -^ (- 1 )^ r(i/2-fc) ^ y 

^27r(r/,r/) ^ fc! r(l -/?/2 -/?A:) V2(r7, r;) J 

g ^ (-1)^ r(-l/2-fc) y 

^ 2 ( 77 ,r/) 0 F^ fc! r(l-/?-/ 3 fc) V2(r/,r/)y 


References 

[ADKS96] S. Albeverio, Yu. L. Daletsky, Yu. G. Kondratiev, and L. Streit. 

Non-Gaussian inhnite-dimensional analysis. J. Fund. Anal, 
138(2);311-350, 1996. 


40 






















[Ben03] 

[BK95] 

[Dal91] 

[dFHSW97] 

[DielO] 

[DODS08] 

[dSElS] 

[EK04] 

[EMOT55] 


C. Bender. An ltd formula for generalized functionals of a 
fractional Brownian motion with arbitrary Hurst parameter. 
Stochastic Process. Appl, 104(1);81-106, 2003. 

Yu. M. Berezansky and Yu. G. Kondratiev. Spectral Methods in 
Infinite-dimensional Analysis. Volume 2, volume 12/2 of Math¬ 
ematical Physics and Applied Mathematics. Kluwer Academic 
Publishers, Dordrecht, 1995. Translated from the 1988 Russian 
original by P. V. Malyshev and D. V. Malyshev and revised by 
the authors. 

Yu. L. Daletsky. Biorthogonal analogue of Hermite polynomials 
and the inversion of the Fourier transform with respect to a 
non-Gaussian measure. Func. Anal. App., 25(2);138-140, 1991. 

M. de Faria, T. Hida, L. Streit, and H. Watanabe. Intersection 
local times as generalized white noise functionals. Acta Appl. 
Math., 46(3):351-362, 1997. 

K. Diethelm. The Analysis of Fractional Differential Equations, 
volume 2004 of Lecture Notes in Mathematics. Springer-Verlag, 
Berlin, 2010. 

G. Drumond, M. J. Oliveira, and J. L. Da Silva. Intersection 
local times of fractional Brownian motions with H G (0,1) 
as generalized white noise functionals. In G. G. Bernido and 
V. G. Bernido, editors. Stochastic and Quantum Dynamics of 
Biomolecular Systems, volume 1021 of AIP Conference Proceed¬ 
ings, pages 34-45, Melville, NY; American Institute of Physics 
(AIP), 2008. AIP Gonference Proceedings 1021. 

J. L. da Silva and M. Erraoui. Generalized grey Brownian mo¬ 
tion local time; Existence and weak approximation. Stochastics, 
87(2);347-361, 2015. 

S. D. Eidelman and A. N. Kochubei. Gauchy problem for frac¬ 
tional diffusion equations. J. Differential Equations, 199(2);211- 
255, 2004. 

A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. 
Higher Transcendental Functions. Volume III. McGraw-Hill 
Book Gompany, Inc., New York-Toronto-London, 1955. 


41 



[Feiri] 

[FKN90] 

[Fox61] 

[GJRdSlS] 

[GKS99] 

[GKSOO] 

[GKU99] 

[GS64] 

[GS68] 

[GV64] 

[HKPS93] 


W. Feller. An Introduction to Probability Theory and Its Ap¬ 
plications. Volume II. Second edition. John Wiley & Sons Inc., 
New York, 1971. 

K.T. Fang, S. Kotz, and K.W. Ng. Symmetric multivariate and 
related distributions, volume 36 of Monographs on Statistics and 
Applied Probability. Ghapman and Hall, Ltd., London, 1990. 

G. Fox. The G and H functions as symmetrical Fourier kernels. 
Trans. Amer. Math. Soc., 98:395-429, 1961. 

M. Grothaus, F. Jahnert, F. Riemann, and J. L. da Silva. 
Mittag-Leffler analysis 1; Gonstruction and characterization. J. 
Fund Anal, 268(7);1876-1903, 2015. 

M. Grothaus, Yu. G. Kondratiev, and L. Streit. Regular gen¬ 
eralized functions in Gaussian analysis. Infin. Dimens. Anal. 
Quantum Probab. Relat. Top., 2(l):l-25, 1999. 

Scaling limits for the solution of Wick type Burgers equation. 
Random Oper. Stoch. Equ., 8(1): 1-26, 2000. 

M. Grothaus, Yu. G. Kondratiev, and G. F. Us. Wick calculus 
for regular generalized stochastic functionals. Random Oper. 
Stoch. Equ., 7(3):263-290, 1999. 

1. M. Gel’fand and G. E. Shilov. Generalized functions. Vol. I: 
Properties and Operations. Translated by E. Saletan. Academic 
Press, New York-London, 1964. 

1. M. Gel’fand and G. E. Shilov. Generalized Functions. Vol. 

2. Spaces of fundamental and generalized functions. Translated 
from the Russian by M. D. Friedman, A. Feinstein and G. P. 
Peltzer. Academic Press, New York-London, 1968. 

1. M. Gel’fand and N. Ya. Vilenkin. Generalized Functions. Vol. 
f: Applications of harmonic analysis. Translated by A. Fein¬ 
stein. Academic Press, New York - London, 1964, 1964. 

T. Hida, H.-H. Kuo, J. Potthoff, and L. Streit. White Noise. 
An Infinite Dimensional Calculus, volume 253 of Mathematics 
and its Applications. Kluwer Academic Publishers Group, Dor¬ 
drecht, 1993. 


42 



[H0UZ96] 

[Ito88] 

[KK99] 

[KLP+96] 

[KLS96] 

[Koc90] 

[KS93] 

[KSS97] 

[KST06] 

[KSWY98] 

[Kuo96] 

[LV14] 


H. Holden, B. 0ksendal, J. Ub0e, and T. Zhang. Stochastic 
partial differential equations. Probability and its Applications. 
Birkhauser Boston Inc., Boston, MA, 1996. 

Y. Ito. Generalized Poisson functionals. Probab. Theory Related 
Fields, 77(l):l-28, 1988. 

N. A. Kachanovsky and S. V. Koshkin. Minimality of Appell-like 
systems and embeddings of test function spaces in a generaliza¬ 
tion of white noise analysis. Methods Funct. Anal. Topology, 
5(3);13-25, 1999. 

Yu. G. Kondratiev, P. Leukert, J. Potthoff, L. Streit, and 
W. Westerkamp. Generalized functionals in Gaussian spaces: 
the characterization theorem revisited. J. Funct. Anal., 
141(2);301-318, 1996. 

Yu. G. Kondratiev, P. Leukert, and L. Streit. Wick calculus in 
Gaussian analysis. Acta Appl. Math., 44(3);269-294, 1996. 

A. N. Kochubei. Diffusion of fractional order. Differential equa¬ 
tions, 26(4):485-492, 1990. 

Yu. G. Kondratiev and L. Streit. Spaces of white noise distri¬ 
butions: constructions, descriptions, applications. I. Rep. Math. 
Rhys., 33(3):341-366, 1993. 

Yu. G. Kondratiev, J. L. Silva, and L. Streit. Generalized Appell 
systems. Methods Funct. Anal. Topology, 3(3):28-61, 1997. 

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo. Theory and 
Applications of Fractional Differential Equations, volume 204 
of North-Holland Mathematics Studies. Elsevier Science B.V., 
Amsterdam, 2006. 

Yu. G. Kondratiev, L. Streit, W. Westerkamp, and J. Yan. Gen¬ 
eralized functions in inhnite-dimensional analysis. Hiroshima 
Math. J., 28(2):213-260, 1998. 

H.-H. Kuo. White Noise Distribution Theory. Probability and 
Stochastics Series. GRG Press, Boca Raton, FL, 1996. 

Joachim Lebovits and Jacques Levy Vehel. White noise-based 
stochastic calculus with respect to multifractional Brownian mo¬ 
tion. Stochastics, 86(1):87-124, 2014. 


43 



[Mai95] 

[Mai96a] 

[Mai96b] 

[Mis08] 

[ML05] 

[MM09] 

[MMPIO] 

[MP08] 

[MS69] 

[MSHIO] 

[OB09] 

[OSSll] 


F. Mainardi. The time fractional diffusion-wave equation. Izv. 
Vyssh. Uchebn. Zaved. Radiofiz., 38(l-2);20-36, 1995. 

F. Mainardi. Fractional relaxation-oscillation and fractional 
diffusion-wave phenomena. Chaos Solitons Fractals, 7(9);1461- 
1477, 1996. 

F. Mainardi. The fundamental solutions for the fractional 
diffusion-wave equation. Appl. Math. Lett., 9(6);23-28, 1996. 

Y. S. Mishura. Stochastic Calculus for Fractional Brownian Mo¬ 
tion and Related Processes, volume 1929 of Lecture Notes in 
Mathematics. Springer-Verlag, Berlin, 2008. 

G. Mittag-Leffler. Sur la reprfeentation analytique d’une 
branche uniforme d’une fonction monogene (cinquieme note). 
Acta Math., 29(1);101-181, 1905. 

A. Mura and F. Mainardi. A class of self-similar stochastic pro¬ 
cesses with stationary increments to model anomalous diffusion 
in physics. Integral Transforms Spec. Funct., 20(3-4);185-198, 
2009. 

F. Mainardi, A. Mura, and G. Pagnini. The M-Wright function 
in time-fractional diffusion processes: a tutorial survey. Int. J. 
Differ. Equ., pages Art. ID 104505, 29, 2010. 

A. Mura and G. Pagnini. Gharacterizations and simulations of 
a class of stochastic processes to model anomalous diffusion. J. 
Rhys. A, 41(28):285003, 22, 2008. 

A. M. Mathai and R. K. Saxena. Distribution of product and 
the structural set up of densities. Ann. Math. Statist., 40:1439- 
1448, 1969. 

A. M. Mathai, R. K. Saxena, and H. J. Haubold. The H- 
Function. Theory and Applications. Springer, New York, 2010. 

E. Orsingher and L. Beghin. Fractional diffusion equations and 
processes with randomly varying time. Ann. Probab., 37(1):206- 
249, 2009. 

M. J. Oliveira, J. L. Silva, and L. Streit. Intersection local times 
of independent fractional Brownian motions as generalized white 
noise functionals. Acta Appl. Math., 113(1):17-39, 2011. 


44 



[PBM90] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev. Integrals 
and Series. Volume 3. Gordon and Breach Science Publishers, 
New York, 1990. More special functions. Translated from the 
Russian by G. G. Gould. 

[Pol48] H. Pollard. The completely monotonic character of the Mittag- 
Leffler function Ea{—x). Bull. Amer. Math. Soe., 54:1115-1116, 
1948. 

[RS72] M. Reed and B. Simon. Methods of Modern Mathematieal 
Physies. I. Funetional Analysis. Academic Press, New York, 
1972. 

[Sch71] H. H. Schaefer. Topologieal Vector Spaces. Springer-Verlag, New 
York, 1971. Third printing corrected. Graduate Texts in Math¬ 
ematics, Volume 3. 

[Sch90] W. R. Schneider. Fractional diffusion. In Dynamics and stochas¬ 
tic processes (Lisbon, 1988), volume 355 of Lecture Notes in 
Phys., pages 276-286. Springer, New York, 1990. 

[Sch92] W. R. Schneider. Grey noise. In Ideas and methods in mathemat¬ 
ical analysis, stochastics, and applications (Oslo, 1988), pages 
261-282. Gambridge Univ. Press, Gambridge, 1992. 

[SKM93] S. G. Samko, A. A. Kilbas, and O. I. Marichev. Fractional In¬ 
tegrals and Derivatives. Gordon and Breach Science Publishers, 
Yverdon, 1993. 

[Sko] A. V. Skorohod. Integration in Hilbert Space, volume 79 of 

Ergebnisse der Mathematik und ihrer Grenzgebiete. 

[ST94] G. Samorodnitsky and M. S. Taqqu. Stable Non-Gaussian Ran¬ 
dom Processes. Stochastic Modeling. Ghapman & Hall, New 
York, 1994. Stochastic Models with Infinite Variance. 

[SW89] W. R. Schneider and W. Wyss. Fractional diffusion and wave 
equations. J. Math. Phys., 30(1);134-144, 1989. 

[Wim05] A. Wiman. Uber den Fundamentalsatz in der Teorie der Funk- 
tionen E'^{x). Acta Math., 29(1);191-201, 1905. 


45 



